o 

< 



00 



o 






A Mathematical Model for the Dynamics and 
Synchronization of Cows 

Jie Sun* Erik M. Bollt^ Mason A. Porter* 



O 
(N 

Ch ' Marian S. Dawkins^ 

^ June 14, 2011 

m 



Abstract 



■^^ ' We formulate a mathematical model for daily activities of a cow (eat- 

ing, lying down, and standing) in terms of a piecewise affine dynamical 
system. We analyze the properties of this bovine dynamical system rep- 
resenting the single animal and develop an exact integrative form as a 
discrete-time mapping. We then couple multiple cow "oscillators" to- 
gether to study synchrony and cooperation in cattle herds. We comment 

Cn ' on the relevant biology and discuss extensions of our model. With this 

^ I abstract approach, we not only investigate equations with interesting dy- 

namics but also develop interesting biological predictions. In particular, 
our model illustrates that it is possible for cows to synchronize less when 
the coupling is increased. 
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1 Introduction 

The study of collective behavior — whether of animals, mechanical systems, or 
simply abstract oscillators — has fascinated a large number of researchers from 
observational zoologists to pure mathematicians |39 1 l47j . In animals, for exam- 
ple, the study of phenomena such as flocking and herding now involves close 
collaboration between biologists, mathematicians, physicists, computer scien- 
tists, and others (TUlIIlllMllil] • This has led to a large number of fundamental 
insights — for example, bacterial colonies exhibit cooperative growth patterns [5] , 
schools of fish can make collective decisions |48], army ants coordinate in the 
construction of bridges [M], intrinsic stochasticity can facilitate coherence in 
insect swarms J52] , human beings coordinate in consensus decision making |18j , 
and more. It has also led to interesting applications, including stabilization 
strategies for collective motion [41] and multi- vehicle flocking [9]. 

Grazing animals such as antelope, cattle, and sheep derive protection from 
predators by living in herds J19U29J . By synchronizing their behavior (i.e., 
by tending to eat and lie down at the same time), it is easier for the ani- 
mals to remain together as a herd [TT1I3D]. When out at pasture, cattle are 
strongly synchronized in their behavior [6], but when housed indoors during 
the winter, increased competition for limited resources can lead to increased 
aggression [Tl[29l[33], interrupted feeding or lying [7], and a breakdown of syn- 
chrony |30| . There is a growing body of evidence that such disruptions to 
synchrony (in particular, disruptions to lying down) can have significant effects 
on cattle production (i.e., growth rate) and cattle welfare P01|21[[^5H271I301I5T| . 
Indeed, synchrony has been proposed as a useful measure of positive welfare in 
cattle [20l|32] , and the European Union regulations stipulate that cattle housed 
in groups should be given sufRcient space so that they can all lie down simul- 
taneously (Council Directive 97/2/EC). In the winter, cattle have to be housed 
indoors; space for both lying and feeding is thus limited, and welfare problems 
can potentially arise because such circumstances interfere with the inherent in- 
dividual oscillations of cows. 

Although cattle synchronize their behavior if space and resources allow, the 
mechanism by which they do this is not fully understood |li p 32 | . In this paper, 
we examine interacting cattle using a mathematical setting to try to gain an 
understanding of possible mechanisms. Viable approaches to studying inter- 
acting cows include agent-based models as well as further abstraction via the 
development and analysis of appropriate dynamical systems to model the cattle 
behavior. In a recent dissertation [22], B. Franz modified the animal behavior 
model of Ref. [l3] to develop an agent-based model of beef cattle and conduct a 
preliminary investigation of its synchronization properties. Given the extreme 
difficulty of actually understanding the mechanisms that produce the observed 
dynamics in such models, we have decided instead to take a more abstract ap- 
proach using dynamical systems. 

Cattle are ruminants, so it is biologically plausible to view them as oscilla- 
tors. They ingest plant food, swallow it and then regurgitate it at some later 
stage, and then chew it again. During the first stage (standing/feeding), they 
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stand up to graze, but they strongly prefer to lie down and 'ruminate' or chew 
the cud for the second stage (lying/ruminating). They thus oscillate between 
two stages. Both stages are necessary for complete digestion, although the du- 
ration of each stage depends on factors such as the nutrient content of the food 
and the metabolic state of the animal [35] |1| We thus suppose that each cow is 
an oscillator, and we choose each oscillator to be a piecewisc affine dynamical 
system in order to incorporate the requisite state-switching behavior in the sim- 
plest possible fashion. Even with this simple model, each individual cow exhibits 
very interesting dynamics, which is unsurprising given the known complexities 
of modeling piecewise smooth dynamical systems [SUTBIEH] . Piecewise smooth 
systems have been employed successfully in numerous applications — especially 
in engineering but occasionally also in other subjects, including biology [231124] . 
To our knowledge, however, this paper presents the first application of piecewise 
smooth dynamical systems to animal behavior. 

Our contributions in this paper include the development of a piecewise affine 
dynamical system model of a cow's eating, lying down, and standing cycles; an 
in-depth analysis of the mathematical properties of this model; investigation of 
synchronization in models (which we call herd models) produced by coupling 
multiple copies of the single cow model in a biologically-motivated manner; 
and a discussion of the biological consequences of our results. Although our 
approach is abstract, the present paper is not merely an investigation of equa- 
tions with interesting dynamics, as we have also developed interesting biological 
predictions. 

The rest of this paper is organized as follows. In Section [21 we discuss the 
dynamical system that we use to describe the behavior of a single cow. We 
present, in turn, the equations of motion, conditions that describe switching 
between different states (eating, lying down, and standing), and a discrete rep- 
resentation using a Poincarc section. In Section [3] we analyze this single cow 
model by studying its equilibrium point, periodic orbits, and bifurcations. We 
examine interacting cows in Section [4] We present the coupling scheme that we 
use to construct our herd equations, introduce the measure of synchrony that 
we employ, and examine herd synchrony numerically first for a pair of cows and 
then for larger networks of cows. In Section [5] we comment on our results and 
bricfiy discuss variant herd models that can be constructed with different types 
of coupling. We then conclude in Section [6] and provide details of our Poincare 
section and map constructions and analysis in Appendix [Al 



^This oscillating approach to eating is one of the things that made cattle suitable for 
domestication, as they can eat during the day and then be locked up safely at night to 
ruminate). 
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2 Single Cow Model 

2.1 Equations of Motion 

We construct a caricature of each cow by separately considering the observable 
state of the cow (eating, lying down, or standing) and its unobservable level of 
hunger or desire to lie down, which can each vary between and 1. We also 
need a mechanism to switch between different states when the level of hunger or 
desire to lie down exceeds some threshold. We therefore model each individual 
cow as a piecewise smooth dynamical system |16j . 
We model the biological status of a single cow by 

w= (a;,y;6') e [0,1] X [0,1] X e. (1) 

The real variables x and y represent, respectively, the extent of desire to eat 
and lie down of the cow, and 

9ee^{£,n,S} (2) 

is a discrete variable that represents the current state of the cow (see the equa- 
tions below for descriptions of the states). Throughout this paper, wc will refer 
to as a symbolic variable or a state variable. One can think of the symbolic 
variable as a switch that triggers different time evolution rules for the other 
two variables x and y. 

Wc model the dynamics of a single cow in different states using 

I ^ ^ —a2X , 
{£) Eating state: < (3) 

[y = (3iy- 

{X ^^^ Oil X 
(4) 
y^ -P2y- 

{S) Standing state: I ' (5) 

where the calligraphic letters inside parentheses indicate the corresponding val- 
ues of 9. For biological reasons, the parameters ai, a2, /3i, and ^2 must all be 
positive real numbers. They can be interpreted as follows: 

ai : rate of increase of hunger , 
a2 : decay rate of hunger , 
Pi : rate of increase of desire to lie down , 
^ /32 : decay rate of desire to lie down . 

The monotocity in each state (growth versus decay) is the salient feature 
of the dynamics, and we choose a linear dependence in each case to facilitate 
analytical treatment. The piecewise smooth dynamical system describing an 
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individual cow is thus a piecewise affine dynamical system [16| . As we shall see in 
the following sections, this simple model is already mathematically interestingo 
Additionally, note that we could have added an additional positive param- 
eter e <C 1 to each equation to prevent the degeneracy of the {x, y) = (0, 0) 
equilibrium point that occurs for all three equationso 

2.2 Switching Conditions 

The dynamics within each state do not fully specify the equations governing a 
single cow. To close the bovine equations, we also need switching conditions 
that determine how the state variable changes. We illustrate these switching 
conditions in Fig. [1] and describe them in terms of equations as follows: 

8 if 61 e {7^, S} and .t = 1 , 

7^ if 61 e {£, 5} and .T < 1 , 2/ = 1 , (6) 

S lie e{£,n] anAx <l,y = 5 {OT x = 5,y<\). 

The positive number 5 < 1 allows the point {x, y) ~ (0, 0) to be excluded from 
the domain, so that the degenerate equilibrium at that point becomes a so- 
called virtual equilibrium, point (i.e., an equilibrium point that is never actually 
reached by the system) [16]. 

Equations (O HI [SJ [6]) form a complete set of equations describing our single 
cow model. This bovine model is a piecewise smooth dynamical system, to 
which some important elements of the traditional theory for smooth dynamical 
systems do not apply, as discussed in depth in the recent book |16) . 

2.3 Discrete Representation 

Although it is straightforward to solve Eqs. (J3l lU [5]) for the fixed state 9, it is 
cumbersome to use such a formula to obtain analytical expressions when the flow 
involves discontinuous changes in 9 (as specified by the switching conditions). 
Therefore, we instead study the dynamics on the boundaries as discrete maps 
rather than the flow on the whole domain. We accomplish this by appropriately 
defining a Poincare section [38] as the surface 

E = {{x,y]9)\x^l,5<y<l,9 = £}\j{{x,y]9)\5<x<l,y = l,9 = 'R} 
= dSDdTZ, (7) 



■^Any differential equation whose flow in a given region (increasing versus decreasing) is 
monotonic in both x and y in aU of the states can be treated similarly using the method we 
describe in Section 2.3 through an appropriate Poincare section. It is expected to produce 
qualitatively similar results, as the detailed flow between state transitions is irrelevant once 
the intersections with Poincare section have been determined. 

^This degeneracy can also be conveniently avoided by restricting the dynamics of x and y 
to a region that excludes the point (0, 0). We opt for the latter choice (see the next subsection 
for details). 
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Figure 1: (Color online) Switching conditions for the single cow model. In the 
left panel, we project the set [S, 1] x [6, 1] x on R^, where edges of the square 
correspond to the borders at which switching occurs. In the right panel, we show 
the detailed switching situations; an arrow from one edge to another indicates 
the change of 9 at that edge from one state to the other. (The arrows with solid 
curves are the ones that leave state 7?,, those with dashed curves are the ones 
that leave state £, and those with dotted curves are the ones that leave state 
S.) 
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which is transverse to the flow of Eqs. ([31 SI [S]) as long as ai^2 and /3i,2 > 0. 
(See the Appendix for the proof.) Furthermore, any flow for which all four of 
these parameters are positive intersects E recurrently (again see the Appendix). 
Although E itself is sufficient to construct a Poincare map (we will use / 
to represent this map on E), it is convenient to consider the discrete dynamics 
on an extended Poincare section E', which we define by adding the other two 
boundaries of the projected square to E to obtain 

E' = i:u{{x,y;e)\x^5,S<y<l}U{{x,y;s)\S<x<l,y = 5} 

= dSUdnu dSy U dS^ , (8) 

where dSx and dSy are used to represent the sets {{x, y;9)\x = S,S < y < 1} and 
{{x,y;9)\S < x < 1 ,y = 5}, respectively. We illustrate the extended Poincare 
section in the left panel of Fig. [H 

The Poincare map on E' is given by the discrete dynamics g : E' -^ E' 
derived by solving Eqs. ([3l |4l [5]) with respect to appropriate initial conditions. 
As we show in the Appendix, this map is given explicitly by 



g{x = l,S<y<l;£) 
g{S<x<l,y = l;n)- 
g{x = S,S <y <1;S) -- 



J {y^' A]T^), if y > S"^ , case (a); 

[{S,S °2y;5), ify<5°2, case(&); 

f i2_ sn_ 

|(l,a;°i;f), ii x > S '^^ ^ case (c) ; 

[{6'T^x,S;S), iix<ST^, case (d) ; 

t>i t>i 

(1,(5 °ij/;£), iiy<5"i, case (e) ; 

-21 il 

. {y ^^ S, 1; 7^) , if y > (5°i , case (/) ; 



giS<x<l,y^S;S) = i^^^ir^'^^^^ if^^^V ^^'^^9); (g) 

[{S ^x, l;7^), iix<S''i, case (h) . 

In Fig. [21 we show all possible mappings on E' and, in particular, illuminate all 
of the possible cases in ^. The Poincare map / : E ^ E can be obtained from 
g (see the discussion in the Appendix). 

3 Analysis of the Single Cow Model 

In this section, we summarize a few properties of the single cow model in terms 
of the discrete dynamics / on E. Specifically, we give analytical results for the 
emergence and stability of the fixed point (which is unique) and the period-two 
orbits on E. We include detailed derivations of these results in the Appendix. 
We summarize these results in Table [1] 

For convenience, we assume that the cow is initially in the state £ with .t = 1 
and 6 < y < 1. If the cow were to start in other situations, it would eventually 
come to this state. Furthermore, we have chosen to assign the state value 9 = S 
to the point (x, y) = (l, 1) for as a tie-breaker. Similarly, 9 = £ ai {x, y) = (1, S) 
and 9 ^ TZ at {x, y) = ((5, 1), in accordance with Eq. (jH]). 
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Figure 2: (Color online) All of the possible rules for determining the discrete 
dynamics on S' that are derived from the original system. For example, from 
9 = £, the flow is going to either hit the horizontal y = 1, which triggers the state 
d ^ TZ [case (a)], or hit the vertical x = S, resulting in the transition 6 -^ S 
[case (6)]. The other three panels similarly demonstrate the other switching 
possibilities for the variable 9. 
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3.1 Fixed Point 

The only possible fixed point on S is the corner point (a;, y; s) ~ (1, 1; £). This 
fixed point is asymptotically stable if and only if the parameters satisfy 

^.f<l. (10) 

Additionally, (when the above condition holds) numerical simulations indicate 
that the basin of attraction of this fixed point seems to be the entire domain. 

3.2 Period-Two Orbits 

The next simplest type of orbits for the discrete map have period two and 
correspond to cycles of the fiow. A period-two orbit on E must contain points 
for which 9 — £ and 9 = TZ appear alternatively. This can occur in a few different 
situations (see Fig. [3|), which we summarize in the following subsections. We 
include further details in the Appendix. Note that some of the period-two 
orbits correspond to higher-period orbits of the discrete dynamics on S'. For 
convenience, we represent such orbits on E', with the understanding that when 
restricted to S (i.e., when points with symbolic variable S are excluded), they 
all have period two. 

3.2.1 Case A: {xo,yo;£) ^ (a■l,yl;7^) -^ {xo,yo;S) ^ ■■■ 
The existence of such an orbit requires the parameters to satisfy 

^4 = 1. (11) 

This implies that there are infinitely many stable but not asymptotically stable 
period-two orbits. The initial value of cc is Xq = 1 and the initial value of y 
(called 2/0, of course) is in the range 

max f (5, (5 °2 = (5 °i 1 < yo < 1 • (12) 

If 2/0 is outside of this range, then one can see using numerical simulations that 
this orbit will necessarily contain a point in this range that becomes a stable 
period-two orbit. 

3.2.2 Case B: {xo,yo;£) -^ {xi,yi;TZ) -^ {x2,y2;Sx) -^ {xo,yo;£) ^ ••■ 
The parameters need to satisfy 

^4>1 (13) 

ai Pi 

or else any trajectory either converges to a fixed point or a stable period-two 
orbit (as discussed above). It is also necessary that 

1111 

— + — > — + — if /3i < a2 . (14) 

ai a2 Pi P2 
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dSj^ (x3,?/3) 



Figure 3: (Color online) Illustration of all of the possible period-two orbits on 
E. 
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There is only one period-two orbit of this type if the above two conditions hold. 
This implies that xq = 1 and 

yo = S'+'-^/o.^ . (15) 

This periodic orbit is asymptotically stable if and only if 

a2 < ai . (16) 

That is, the orbit is asymptotically stable if and only if the rate at which a 
cow becomes sated while it eating is slower than the rate at which it becomes 
hungrier when it is not eating. 

3.2.3 Case C: {xo,yo;£) -^ {xi,yi;Sy) -> {x2,y2;T^) -^ {xo,yo;£) -^ ■ ■ ■ 
Again, we first need 

ai Pi 
Additionally, 

— + —<— + — and Pi <a2. (18) 

ai a2 Pi P2 

There is also only one period-two orbit of this type; it has xq ~ 1 and 

l/ai + l/a2 
yo := 51/-91+1/-92 . (19) 

This orbit is asymptotically stable if and only if 

I32<l3i. (20) 

This case is analogous to case B, except that the roles of lying down and eating 
have been reversed. Hence, this period-two orbit is asymptotically stable if and 
only if the rate at which a cow desires to get up when it is lying down is slower 
than the rate at which it increases its desire to lie down when it is not lying 
down. 

3.2.4 Case D: {xo,yo':£) -> {xi,yi;Sy) -^ {x2,y2]Tl) -^ (3:3,^3; 5^;) -J- 
(2:0,2/0;^) ^ ■■■ 

The appearance of this orbit requires the following conditions to be satisfied: 

'^_h > 1 
ai Pi 

— + — = — + — and /3i < a2 . (21) 

ai a2 pi P2 

There are infinitely many such orbits, which satisfy xq = 1 and 

S <yo<S^ . (22) 

All of these orbits are stable but not asymptotically stable. 
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Table 1: Summary of low-period orbits (up to period two) and their stability of 
the single cow dynamics restricted to the Poincare section E. All orbits except 
for the first one are period-two orbits on E. In the 'Stability' column, we use 
'a.s' as an abbreviation for 'asymptotically stable'. 



Parameters 


Orbit 


Condition on j/o 


Stability 


^.^<1 
ai Pi 


{(1,1;^)} 


none 


a.s 


a2 P2 _^ 
ai Pi 


{{l,yo:£),{yJ^,l-n)} 


/3l 

max ((5, (5 °2 ) < yQ < 1 


stable 


— ■ — >l,a2< Pi 
ai Pi 


{(i,yo;^),(#,l;7e)} 




a.s iff a2 < Oil 


I ai a2 ~ Pi P2 


{(i,yo;^),(#,i;7e)} 


2/0 = d °i 


a.s iff a2 < Oil 


— • — > 1,^2 > /3i; 
I ai Pi 

S 1 1 1 1 

lai a2 Pi p2 


{(l,yo;f),(<5,(5-^^/o;7^)} 




a.s iff P2 < Pi 


— • — > 1,02 > /3i; 
I ai Pi 

S 1 1 11 

I ai a2 Pi P2 


{(l,yo;f),(<5'+-yo~^>l;7^)} 


(5 < j/o < ^ °2 


stable 



3.2.5 Summary 

Wc summarize the emergence of low-period orbits (up to period two) of / : E — >■ 
E in different parameter ranges in Table [TJ 

3.3 Grazing Bifurcations 

We remark that the single cow equations cannot exhibit grazing bifurcationso 



3.4 Higher-Period Orbits and Bifurcation Diagram 

Although one could proceed to analyze more complicated orbits, this is not the 
main topic of this paper. Instead, we simply illustrate the existence of more 
complicated (possibly chaotic) orbits through a bifurcation diagram (see Fig. [J]) 
by simulation with varying one of the parameters. This parameter, which we 
choose to be a2 , seems to be transverse to the unfolding of the bifurcation and 
reveals rich dynamics in our model. 

For a wide range of parameters, there seems to always be a dense subset (for 
a fixed set of parameters) of the domain that attracts "typical" (in the sense 

*In the theory of piecewise smooth dynamical systems, a grazing bifurcation is said to 
occur when a limit cycle of a flow becomes tangent to a discontinuity boundary [8l ll6| . 
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Figure 4: (Color online) Bifurcation diagram for the discrete dynamics / on S. 
We fix the parameter values ai = 0.05, /3i = 0.05, /32 = 0.125, and i5 = 0.25. 
The vertical axis is ^ = j/ + (1 — x), corresponding to the points {x,y) on S. 
In the top panel, we show the diagram for which q = ^^ ranges from to 5; 
dashed and dotted curves give theoretical results, which we summarize in Table 
[TJ In the bottom panel, we show the diagram for q from to 15. If we further 
increase q, the two large finger-like bands on the right of the diagram retain their 
shape and become progressively closer. Numerical simulations suggest that the 
distance between them tends to as gr — > oo. 
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of nonzero measure) initial conditions. We show one of these (hkely chaotic) 
orbits in Fig. [5l We connect the dots by straight lines in order to illustrate the 
end points of the flow touching the boundaries, although the actual trajectories 
between points on the boundaries are convex curves. We remark that one can 
think of the discrete dynamics on E' as a billiard-like problem (see Refs. [T5ll45] 
and references therein for discussions of billiards) with nontrivial bouncing rules 
on the boundary and nonlinear potentials that determine the trajectories of 
particles between collisions with the boundary. 



0.75 




0.25 



Figure 5: (Color online) A typical discrete orbit (thin solid lines) on S for 
the parameters ai = 0.05, ^2 = 0.1, /3i = 0.05, ^2 = 0.125, and S = 0.25. 
We depict the case corresponding to g = 2 in Fig. U] The dashed lines show 
transient dynamics. We highlight the boundaries using thick solid lines. For 
aesthetic reasons, we join successive points on S with straight lines, and we note 
that the actual flow that connects these points are piecewise convex curves. 



4 Coupled Cows and Synchronization 

As we discussed in the introduction, there are many biological benefits to achiev- 
ing synchronized eating and lying down in cattle. We are thus motivated to con- 
struct herd equations that describe interacting cows by coupling the single cow 
equations (jSHHl)- We make specific choices motivated by biology and simplicity, 
though it is of course important to consider both more complicated choices and 
alternative forms of coupling. Our goal is to highlight just one possible form of 
the interactions in detail, but we hope that our work will serve as a springboard 
for rumination of some of the alternatives that we will mention briefiy in Section 

m 
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In this section, we numerically investigate the effect of coupling in a system 
of a few cows. For the purpose of simplifying the exposition of the equations, 
we use indicator functions defined on the set {£,TZ,S}: 

Xi,iO) = < . (23) 

I , otherwise . 

The single cow equation in between state transitions ^ can then be written 
as 

{x = a{6)x, (24) 

where we have defined functions 

\ot{9) = -a2X£{0) + oiixnid) + aixs{0) , /25) 

\m ^ l^ixm ~ hxn{0) + PiXs{e) . ^ ' 

4.1 Coupling Scheme 

There arc numerous possible ways to model the coupling between cows. Wc have 
chosen one based on the hypothesis that a cow feels hungrier when it notices the 
other cows eating and feels a greater desire to lie down when it notices other 
cows lying down. (Wc briefly discuss other possibilities in Section [5l) This 
provides a coupling that does not have a spatial component, in contrast to the 
agent-based approach of Ref . [55] . We therefore assume implicitly that space is 
unlimited, so we are considering cows to be in a field rather than in a pen. We 
suppose that the herd consists of n cows and use i to represent the z-th cow in 
the herd. This yields herd equations given by 

{x, ^ [aW(sj) + ^Y.%l'^^jXe{sj)\x, , 

ly.= [/3«(s.) + tE;=ia,,X7^(s,)]y,, ^ ^ 

with switching condition according to Eq. ([6|) for each individual cow. The 
second terms in both equations give the coupling terms of this system. The 
matrix A = [aij\nyin is a time-dependent adjacency matrix that represents the 
network of cows. Its components are given by 

I 1 if the i-tli cow interacts with the j-th cow at time t , 
aij [t) = \ , 

I if the i-tli cow does not interact with the j'-th cow at time t . 

(27) 
Additionally, ki = X^jLi ^ij is the degree of node i (i.e., the number of cows to 
which it is connected), and the coupling strengths ax and ay are non- negative 
(and usually positive) real numbers corresponding to the strength of coupling. 
This is designed to emphasize that animal interaction strengths consider prox- 
imity to neighboring animals. 
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It is important to note that in the case where A is time-independent, the 
dynamics governing the network of interacting cows only changes when at least 
one of the individual cows changes its state 9i. In practice, we can solve ana- 
lytically for the flows in between such transitions (because they are piecewise 
affine differential equations) instead of performing numerical integration in the 
whole time interval, which might cause numerical instability when the number 
of transitions becomes large. 

4.2 Measuring Synchrony 

Wc also need a measure for the synchrony between cows. For each cow i, let 
t'-*' and K^'^ be vectors such that 

{r. = The k-th time at which the i-th cow switches its state to £ , ,„„, 

(i) (28) 

K^, ' = The k-th time at which the i-th cow switches its state to TZ . 

Given pairs of vectors r*-*' and t^^' of the same length, the "eating" synchrony 
between cows i and j is measured by 

K 



(k«--<^'l>=^El-^'--^''l' (29) 



^J M 1/ ^ 

fe=l 



where (•) denotes time-averaging. In general, the vectors t^^' and r'^' are of 
different lengths, so wc truncate and shift one of them to match up with the 
other in such a way that it gives approximately the minimal Af- as defined 
above. 

Similarly, we define the "lying" synchrony between cows i and j by 

A5^(|^«-«:«|). (30) 

For n cows, the group "eating" and "lying" synchrony are then measured by 
averaging over all of the synchrony between individual pairs: 

I A^ EE (Af^.) = :^ J2i,j Afj' , . 

IA^-(A-)=;^E,,A5> ^ ' 

and the aggregate synchrony can then be measured via 

A EE A^ -f A'^. (32) 

There are, of course, other possible measures of synchrony that one could em- 
ploy. For example, in his agent-based study, Franz [52] considered kappa statis- 
tics, an order parameter adapted from the usual one used in the Kuramoto 
model, and a direct count of how often all cows are lying down ^01155] . 

4.3 Numerical Exploration of Herd Synchrony 

With the tools described above, we are now ready to show some examples of 
synchronization of cows. We will start with a system consisting of only two 
cows and then consider herds with more than two cows. 
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4.3.1 Two Coupled Cows 

We first examine how the coupUng strength affects the extent of synchronization. 
Assume that the two cows have individual dynamics that arc specified by nearly 
identical parameters: 



a 



(1,2) 



P\ 



(1,2) 



= 0.05 ±e, 
= 0.05 ±e, 
= 0.25. 



a. 



(1,2) 
2 

,(1,2) 



= 0.1±e, 
= 0.125±e. 



(33) 



We show simulation results in Fig. [5] and Fig. [7] to illustrate the dependence 
of synchrony both on the parameter mismatch e and on the coupling strength 
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Figure 6: (Color online) Typical time series of the state variables 61,2 for different 
coupling strengths. The system of equations is described by Eq. ([26]) . and the 
parameter values are given in Eq. p3p . The parameter mismatch between the 
two cows is e = 10"'^. The horizontal axis is time t. The left panel shows 
the transition of states 61 (red circles connected by '— ') and 62 (black crosses 

connected by ' ') of a typical time series with the coupling strengths a.^. = 

(Ty = (i.e., when there is no coupling). The right panel shows a similar plot 
with the coupling strengths a^ = cFy = 0.045. 

These pictures suggests that our measure of synchrony is reasonable for 
such a system. The greater the difference between the two cows, the harder it 
is for them to achieve synchrony. However, the dependence of synchrony is not 
necessarily monotonically dependent on the coupling strength. An increase in 
the coupling strength at the beginning does improve synchrony, but there is a 
point beyond which larger coupling can in fact lead to lower synchrony. 



4.3.2 Network of Coupled Cows 

In this subsection, we show numerical results on synchronization among a few 
cows. In all examples, we consider a herd of n = 10 cows, where each individual 
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Figure 7: (Color online) Dependence of synchrony on coupling strength. The 
system of equations and parameters are specifed by Eq. ([26| and Eq. ([33|) . 
respectively. In the left panel, we illustrate the synchronization error, which 
we measure using Eq. (PT|) . for different coupling strengths ax,y for two coupled 
cows whose parameter mismatch is e = 10~^. In the right panel, we show the 
synchronization error for parameter mismatch e = 10~^. We obtain each curve 
(for a fixed e) by averaging over 50 runs. Each run is an independent realization 
of the herd equations starting from an initial condition chosen uniformly at 
random. Vertical error bars indicate one standard deviation from the mean. 



has parameter values slightly perturbed from ai — 0.05, a2 = 0.1, /3i — 0.05, 
and (^2 = 0.125. Additionally, we note that 10"'^ is the maximum difference in 
each parameter value relative to the average parameter among all individuals. 
We can couple these cows using different network architectures — for example, 
a circular lattice and a star graph (see Fig. [8]). We use these networks only as 
illustrative examples, as one can of course perform similar investigations with 
any other network architecture. 

In Fig. [9l we show the state transitions of the ten cows during a small time 
interval. We consider fixed coupling strengths cr^ = ctj, = 0.05 for each of the 
two network architectures. 

In Fig. I10[ we illustrate the dependence of synchrony on different coupling 
strengths for the two network configurations. Interestingly, when the coupling 
strength is increased, the cows tend to synchronize less when they are coupled 
via a circular lattice, whereas synchrony is improved if they are coupled via 
a star graph. We have also tested numerically other network configurations, 
such as circular lattices with more than just nearest-neighbor connections and 
(Erdos-Renyi) random graphs, and the resulting curves are qualitatively similar 
to the one shown in the left panel of Fig. [TOl It would be interesting to study 
what network architectures can lead to good synchrony beyond the star graph, 
which is an idealized example. A heuristic reason that synchrony can decrease 
when coupling is increased in our herd model (and, more generally, in piecewise 
smooth dynamical systems) is that decreasing the difference in the observable 
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Figure 8: (Color online) Example network architectures for coupled cows: (left) 
Circular lattice with 10 nodes and (right) star graph with 10 nodes. [The 
spherical cow image was created for this paper by Yulian Ng and used with her 
permission.] 
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Figure 9: (Color online) Typical state transitions for coupled cows in (left) a 
circular lattice and (right) a star graph with fixed coupling strengths a^ ~ Uy = 
0.05. We plot (artificial) straight lines to help visualize transitions between 
states (which are represented by open circles, with different colors representing 
different cows). The horizontal axis is time. Some of the curves overlap (so that 
fewer than 10 colors are visible) due to the partial synchrony between individual 
cows. 
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variables {x, y) through couphng does not necessarily reduce the difference in 
the hidden state variable 9, and the effect of coupling might well be the opposite 
from what one would naively anticipate, as we have observed using the circu- 
lar lattice structure. Moreover, recent work in other contexts has illustrated 
that increasing the number of connections in a network can sometimes lead to 
less synchrony [34| . Although synchronization has been studied extensively for 
smooth dynamical systems [2lgll34 l l37 l [39 l l42 l l46 l l47 l l49 l [50] and the mechanisms 
that promote synchrony in such situations are relatively well understood, lit- 
tle is known about networks of coupled piecewise smooth dynamical systems. 
It would be interesting to study the influence of network architecture on syn- 
chrony for models other than smooth dynamical systems (such as the herd model 
considered in this paper), which might prove to be important in studying the 
behavior of interacting animals. 





0.01 0.02 0.03 0.04 0.05 

a 
x.y 



-0.2L 



0.01 0.02 0.03 0.04 0.05 

o 
x.y 



Figure 10: (Color online) Synchrony measure versus coupling strengths in the 
(left) circular lattice and (right) star graph. 



5 Discussion 



We have only scratched the surface concerning the modeling of herd synchrony 
in cattle. 

We considered each cow as an oscillator, which we modeled as a piecewise 
affine dynamical system. Our single cow model had interesting mathematical 
properties, which we discussed in detail. Monotonic dynamics within each state 
was the most important detail, and we chose affine monotonic dynamics to make 
the analysis as tractable as possible. 

We illustrated herd dynamics through specific coupling choices between 
cows. We assumed that the herd is in a field rather than a pen and, in particular, 
ignored the presence of spatial constraints. We considered cows that become 
hungrier when they notice others eating and a greater desire to lie down when 
they notice others lying down, but numerous other choices would also be inter- 
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esting to study. For example, the relative importances of the two aforementioned 
types of positive coupling can be varied systematically, and the specific func- 
tional forms of coupling can also, of course, be different. Additionally, it is not 
clear whether cow synchrony arises from such an active mechanism or whether 
it can arise from more passive forms of coupling. In particular, this could entail 
the incorporation of spatial effects, such as limited eating and bedding areas 
and the competition of cows for such resources. The inherent oscillations of 
individual cows can lead to synchronization even with almost no interactions 
between individuals [44]. Synchrony can potentially emerge even if the only 
interaction between cows occurs when one steps on another one, so such a mini- 
malist but biologically meaningful mechanism (in which the cows need not even 
notice whether another cow is feeding or resting) would be interesting to test 
against more complicated forms of interaction. It would also be interesting 
to use real observations of cattle to compare the synchronization properties in 
limited space versus "unlimited" space (i.e., pens versus open fields), and such 
experiments are currently in progress. 

One could examine spatial eff^ects in the oscillator model of cows by con- 
sidering more realistic network architectures. Such networks could either come 
from experimental data (which has not yet been gathered) of which cows come 
into contact with each other or using structures that respect the fact that fields 
and pens are planar regions. It would also be interesting to consider different 
network structures from an abstract perspective in order to test observations 
such as the different dynamics with the star graph (which has one high-degree 
node and many small-degree nodes), and also to consider the synchronization 
dynamics of larger herds. Additionally, herds of cattle are known to have hi- 
erarchies, as not all cows are created equal, and this can be incorporated into 
the model either through an appropriate network architecture or by considering 
heterogeneity in the dynamics of individual cows. 

An alternative modeling choice would be to consider agent-based models for 
the herd [22] rather than the oscillator model that we have studied. Agent-based 
formulations are good at incorporating spatial effects, but they of course have 
a black-box fiavor that makes them very difficult to analyze. 

The inherent oscillation between the standing/eating phase and the lying/ruminating 
phase has interesting biological consequences. For example, to stay together as 
a herd, it is not necessary for all cows to be exactly synchronized, as is some- 
times believed. It is possible (and it has been observed often in fields) for a 
herd to have some individuals lying down and other individuals standing and 
grazing around them. From a functional perspective, it is conceivable that this 
could lead to better spotting of predators than if everyone had their heads down 
at the same time. A degree of desychronization (provided that it didn't lead 
to the herd breaking up) might actually be better for each individual than per- 
fect synchronization [3] . Intriguingly the recent model of groups of animals by 
Dostalkova and Spinka, in which each individual can either move or stay in 
one place, has found evidence (using optimization of a cost function) of partial 
synchronization but that completely synchronized and completely desynchro- 
nized situations seem to occur for a much larger set of parameter values |17j . 
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Moreover, their "paradoxical" prediction that average group size might decrease 
as the ratio of the grouping benefit to grouping cost increases is in some sense 
similar (at least philosophically) to our prediction that less synchronization can 
potentially occur even with stronger coupling between individual cows. 

Although we have framed our discussion in terms of cows, our oscillator 
framework is very general and should also be useful — perhaps with modifica- 
tions that are tailored to different species — in studying the behavior of other 
ruminants. It is of considerable biological interest to establish empirically which 
mechanisms for synchrony actually operate in real cows (and, more generally, in 
other ruminants and in other animals) and to discern more precisely the extent 
to which such synchrony actually occurs. It is thus important to develop testable 
predictions that can help one distinguish the numerous possible synchronization 
mechanisms. We have taken one small step in this paper, but there is clearly 
a lot more interesting research on the horizon. It is also desirable to consider 
practical situations, such as the effect of changing pen shape, stocking density, 
size of lying area, feed-trough size and position, and the nutrient quality of the 
food. 

In addition to the many fascinating animal-behavior questions, the research 
reported in this paper also suggests several interesting abstract questions. For 
example, although the theory of synchronization is well-developed and widely 
used for smooth dynamical systems [2l|4l[37l|49l[50], it is an open problem to 
predict in general when a system that is composed of coupled piecewise smooth 
oscillators can achieve a stable synchronous state. In pursuing such considera- 
tions, it would also be relevant to consider different notions of synchrony. Such 
analysis is of potential importance given the wealth of piecewise smooth dy- 
namical systems that arise in many applications J16j . Furthermore, the effects 
of delay and changes in the network architecture in time are also expected to 
affect the synchronization properties, though such considerations arc difficult 
even for smooth systems [42ll43ll46] . We hope that that the model that we have 
developed in this paper will stimulate research along these lines. 

6 Conclusions 

Wc modeled the eating, lying, and standing dynamics of a cow using a piece- 
wise affine dynamical system. We constructed Poincare maps to examine the 
system's equilibrium point and low-period cycles in depth and illustrated more 
complicated behavior using bifurcation diagrams. We then considered a model 
of coupled cows — first using two cows and then using networks of interacting 
cows — in order to study herd synchrony. We chose a form of coupling based 
on cows having an increased desire to eat if they notice another cow eating 
and an increased desire to lie down if they notice another cow lying down. We 
constructed a measure of synchrony that keeps track of when each cow is in 
a given state and showed that it is possible for cows to synchronize less when 
the coupling is increased. We also discussed other forms of coupling and cow- 
interaction networks that can be studied using our formulation. This line of 
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inquiry seems very promising and that it will not only lead to interesting future 
theoretical investigations but can even motivate new experiments. Although we 
framed our discussion in terms of cows, our framework is general and it should 
be fruitful in the study of the behavior of other ruminants as well. The stakes 
are high when studying animal behavior, and we believe that our model of cattle 
herds (and generalizations of our model) will yield increased understanding of 
their synchronization properties. Milking these ideas as much as possible should 
prove to be very insightful from both theoretical and practical perspectives. 
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A Investigation of the Single Cow Model using 
Poincare Section 

The single cow model for w = {x,y;9), with {x,y) S [0,1] x [0,1] and 9 G 
{£,TZ,S}, consists of equations describing dynamics for different states and 
rules for how to switch states. The equations within each state are 

{X ^ — Q^9X 
(34) 
y = Piy, 

(TZ) Resting state: < ' (35) 

[y = -/32y, 

(S) Standing state: < (36) 

[y = /3i2/, 

The rules for switching the state 9 are 

(S if 6'e {7e,5} anda; = 1, 
9^ <n if 6* e {f, 5} and X < 1 , y = 1 , (37) 

[s a 9 e {£,11} a.nd X < I ,y ^ S {or X ^ S ,y < I) . 

All of the parameters (ai,2 and /3i_2) are positive. We use the term single cow 
equations to refer collectively to Eqs. P4l35l36l37p . 

A.l Transversality of the Poincare Section 

As with smooth systems, the term "flow" in piecewise smooth dynamical systems 
designates the usual time-parameterized continuous group |16| . 

Definition 1 (Flow) The solution to the single cow equations, which we de- 
note by (f>{t — to,W()) for initial condition wq at time to, is called a flow of the 
single cow equations. 

The two strips of the boundary of the single cow equations form a set that 
we denote by S. It is defined by 

E = {ix,y;9)\x = l,6<y<l,s = £}U{ix,y;9)\S<x<l,y = l,9 = n} 
= d£Udn, (38) 

where we recall that d£ and dTZ are used to represent the two sets {{x, y; 9)\x = 
1,S <y<l,s^£} and {{x,y;9)\S < x < 1 ,y = 1 ,9 = TZ}. 

The following lemma shows that the surface S can be used as a Poincare 
section for any fiow. This result follows directly from the equations of motion. 

Lemma 1 (Transversality and Recurrence of E) For any initial condition 
Wq = (xo,yo;^o) with initial time to, the flow (f>{t — tQ,WQ) of the single cow equa- 
tions is transverse to T,. In other words, the direction of the flow (restricted to 
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the xy-plane) is not tangent to E (also restricted to the xy-plane). Furthermore, 
there exists t > to such that (j){t — tg, wq) G E. 

A similar lemma holds for the extended Poincarc section E'. which is defined 
as 

E' = j:u{{x,y;0)\x^S,d<y<l}U{{x,y;e)\S<x<l,y^S} 

= d£UdTZU dSy U dS^ , (39) 

where dSx and Sy are used to represent the sets {{x, y;9)\x = S,S < y < 1} and 
{{x,y;9)\S < x < 1 ,y = S}, respectively. 

A. 2 Discrete Dynamics on the Poincare Section: Deriva- 
tion 

The derivation of map g on E' involves first solving for the flows on the contin- 
uous segments where 9 takes one value. 
Starting from 9 — 8, we get 

[ten - ^log(i), g£n{x,y;£) = (2/^,l;7^), 
{x,y;£)^{ ,^ (40) 

[t£s = ^log(i) , gesi^.V.S) - (-5, {\)-y\S) . 

Starting from 9 = TZ, we get 

{ ,82 
tns^ -^^og{^), 57^£(a;,2/;7^) = (l,a;°i;f), 
(41) 
iK5 = ^log(i), gnsix,y;n) = {{^)^2x,d;S). 

Starting from 9 ~ S, we get 

rt5£ -^log(i), gse{x,y;S)^{l,{^)^y-8), 
{x,y;S)~^< ^^ (42) 

U5K = ^log(i), gsn{x,y;S)^{{^)^x,l;TZ). 

Subscripts such as STZ indicate the transition of 9 from one state (e.g., £) 
to another (e.g., TZ). The quantity t with the appropriate subscript represents 
the time it takes for this transition to happen. In the next subsection of this 
appendix, we will analyze the dependence of the discrete system specified by 
the above rules on the parameter values and initial conditions. 

Using the above equations, we derive the discrete dynamics g on E' from the 
nth transition to the (n-|- l)t/i transition. This is given by (x„+i,y„+i, ^n+i) = 



SYNCHRONIZATION OF COWS 26 



g{x„,yn,0„), where 

g{x^l,d<y<l;£) 

g{5<x<l,y^l;TZ)- 

gix^S,S<y<l;S)-^ 



\{y^,l;TZ), iiy>S^ 

-Bx ix 

{{5,5 ^^y;S), ify<(5°2 

'(l,x^;£), \ix>5^ 

^{5^T^x,5]S), ]ix<5^ 

'{l,5~^y-£), ify<(5^ 

Xy-^^5,l-n), ifyxS^ 



^ „i ' ^ ~ »i ' (43) 

((5-'^x,l;7^), iix<5^^. 

This, in turn, yields the discrete dynamics on E'. The dynamics / on S is then 
simply g restricted to S. 

Wc need the following definitions in order to discuss of stability of orbits on 
the dynamics on E. Wc start by defining an appropriate distance measure on 
E. 

Definition 2 (Distance measure on E) We define the distance || • || on E hy 

\\wi-w2\\ = \xi~X2\ + \yi-y2\, (44) 

where Wi = {Xi,yi]9i) for i = 1,2. Note that the symbolic variable 9 does not 
affect the distance. 

We now give a definition of stability and asymptotic stability that is analo- 
gous to the standard definition for smooth dynamical systems [TBI 155] . 

Definition 3 (Stability and Asymptotic Stability) Let f denote the dis- 
crete dynamics on E. A fixed point wq on E is stable if for any e > 0, there 
exists an rj > such that 

\\w-wa\\ <?7^ ||/(«^)-/(wo)ll < £• (45) 

A fixed point wq is asymptotically stable if it is stable and there exists an rj > 
such that 

\\w ~ wo\\ < 7] ^ lim /"(w) = Wo ■ (46) 

n— »oo 

The stability and asymptotic stability of a period-T orbit {zq, . . . , zj^-i} is de- 
fined as the stability and asymptotic stability of the fixed point Zq of the T-th 
iterate f'^ of the map f on Y,. 

A. 3 Fixed Points: Existence and Stability 

We first show that there is only one fixed point on E. 
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Lemma 2 (Fixed Points on S) The only fixed point of f on Yl is the point 
wq — {xo,yo;9) ^ (1,1;£). This fixed point is (locally) asymptotically stable if 
and only if 

(47) 



02 ^ 



Proof 1 First we show that if {xo,yo',Oo) = {l,yo;£), where S < yo < 1 or 
(xcyoS^o) = ixo,l;TZ), then it is not a fixed point. Suppose that there is a 
fixed point starting from wq ~ {xq ~ 1,S < yo < l;do — £). Because it is 
a fixed point on E, the flow cannot hit dTZ. It must thus intersect dSy first 
and then continue and hit d£ again; see Fig. \11\ for an illustration. However, 
because the y-component increases exponentially with rate /3i > when both 
6 = £ and = S — see Eqs. I{34]35138\) — it follows that j/i > j/o- Consequently, 
(xijj/i) cannot be the same point as {xo,yo)- A similar argument applies to 
initial conditions with 9q = TZ, so we can conclude that there is no fixed point 
for the discrete dynamics on S — {(1, 1,£)}. 




Figure 11: (Color online) Illustration that there cannot be any fixed point of 
/ on E except for the corner point (1, 1]£). This follows from the monotonic 
increase of the y-component when 9 = £ remains unchanged. The other pos- 
sible situation (not pictured) occurs when 9 = TZ, for which the x-component 
increases monotonically, indicating that there cannot be an equilibrium point 
on dTZ. 



The only possible fixed point on E is the point (1,1; £). The asymptotic 
stability of this fixed point is easily obtained through linearization. 

We remark that although linearization gives local asymptotic stability of the 
fixed point, numerical simulation indicates that the actual basin of attraction is 
the entire domain when Eq. (|Tfl) is satisfied. 
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A. 4 Period-Two Orbits: Existence and Stability 

We next analyze all possible period-two orbits of / on S. Some of those orbits 
correspond to higher-period orbits of g on E'. When this is the case, we list the 
points of such an orbit on E' to differentiate between different periodic orbits. 
Nevertheless, it is useful to keep in mind that when restricted to E (i.e., when 
one ignores points such that 9 = S), such orbits have period two. Figure [T2l 
illustrates all of the possible period-two orbits. 



{xy,yi)dn 




dS^, ix3,y3) 



Figure 12: (Color online) Illustration of all of the possible period-two orbits on 
E. 



A. 4.1 Case A: {xo,yQ;£) -^ (xl,yl;7^) -^ {xo,yo;£) ^ ■•■ 
This period-two orbit satisfies 

(xo,2/o;^) = (l,2/o;f), where < yo < 1; 



Wo 



Wo 



fi^o) = g{wo) = wi 
P{wo) ^ g^{wo) = W2 



(xl,J/l;7^) = {y^^ ,l;7e), where yo > S' 



fi 



f>2 



(X2,y2]£) = (l,a;i°' ,£), where xi > (5%8) 
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The existence of this orbit entails that (xo,yo) = (3^2,2/2) and that the con- 
straints (i.e., the inequahties that accompany the equations) are satisfied in 
(|i5|) . It is thus required that the parameters satisfy 

^.f = l, (49) 

and that yo satisfy 

(S<yo<l, ifa2</3i, 

< ffi_ (5U) 

[o"^ < yo <l, if a2 > /3i . 

Linearization shows hnear stability of the orbit but not necessarily asympotic 
stability. 

As we are taking xq = 1 (the initial point is on the right edge of the square 
domain), we obtain conditions for yo- All orbits must hit the right edge at some 
point, so we do not lose any generality by taking a;o = 1. 

A.4.2 Case B: {xa,ya;£) -^ {xi,yi;n) -^ {x2,y2]S) -^ {xa,ya;S) ^ ... 
In this case, the orbit has period two on S but period three on S'. Specifically, 

Wo = (-^Ojyoi^o) = (l,yo;'?), where (5 < j/o < 1; 
f{wo) ^ 9{wo) = wi = (.Ti,2/i;6'i) := (ypi ,l;7^), where yo > '^"^ ; 

g'^{wo)^W2 = (.T2,2/2;^2) = ('5~'^xi,(5;5) , where cci < f^*^ ; 
wo = /^(wo) =5^(wo) = W3 = {x3, y3; 63) ^ {1,X2'''S;£), where X2>ST^1) 
For fixed parameter values, there is only one such orbit; it must satisfy 



1+ 



yo = 6'^^. (52) 

The existence of this period-two orbit also requires that the parameters 
satisfy 



(53) 



(54) 

In particular, it is worth remarking that the orbit is not asymptotically stable 
in the case ai — a2 describing equal growth and decay rates for hunger. 
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A. 4. 3 Case C: {xo,yo;S) -^ {xi,yi;S) -^ {x2,y2]Ti) -^ {xo,yo;S) ^ ... 
In this case, the orbit has period two on S but period four on E'. Specifically, 

Wo = {xq, yo; So) = {l,yo;£), where S <yo <1; 

ffi 31 

g{wo) = iui = {xi,yi;si) = {6,6 "2 y^; S) , where yo < 6 •'■2 ; 
f{wo) = g'^iwo) = W2 = {X2,y2;s2) = {y 1^6, l;n), where yi > 6~ ; 

£2. ^ 

Wo = f'^ {wo) = 9^ {wo) = W3 = {x3, y3; S3) = {1,X2^;£), where X2> 61^2 (^55) 



Solving (|55l) with the associated constraints yields necessary conditions for 
the existence of this period- two orbit. The initial value j/o must satisfy 



1 , 1 

— ^^ 



y„^5^^^ , (56) 

and the parameters must satisfy 

02. ^2 -^ -I 

ai ' /?l -^ ' 

ai a2 ^1 & ' 

02. ^ ^ 
[^ ai — I3i ' 

This orbit is asymptotically stable if and only if 

^<1. (58) 

Pi 

Note, in particular, that this implies that the orbit is not asymptotically stable 
when /3i =132 (i.e., when the growth and decay rates for desire to lie down are 
equal) . 

A. 4. 4 CaseD: {xo,yo;£) -^ {xi,yi;S) -> {x2,y2;T^) -^ (2^3, 2/3; 5) -^ {xo,yo;£) 



In this case, the orbit has period two on E but period four on E'. Specifically, 

Wo == (a;o,2/o;so) = (l,yo;^), where (5 < 2/0 < 1; 
g{wo) = wi = (.Ti,2/i;si) = ((5, (5 °-2 y^-^ S) , where yo < 5 "2 -^ 

f (wo) = 9 {wo) = W2 = (x2, 2/2; S2) = (^1 "(5,1; 7^), where 2/1 > (5°i ; 
g^{wo) = W3 = {x3,y3;s3) = {6~^X2,6;S) , where X2 < 6^ ; 

wo = f{wo)= g^{wo) = w^ = (a;4,2/4;s4) = (l,a;3 °i(5;f), whereas > (5^ . 

(59) 
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The existence of such orbits entails that 

ai Pi 
1111 

ai a2 /3i /32 

Pi < a2. (60) 

This yields infinitely many such orbits, for which xq = 1 and 

S<yo<S^ . (61) 

All of these orbits are stable but not asymptotically stable. 
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